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Abstract 



^ . We study the entanglement effect of beam splitter on the temporally stable phase states. Specif- 

I ically, we consider the eigenstatcs (phase states) of an unitary phase operator resulting from the 

polar decomposition of ladder operators of generalized Weyl-Heisenberg algebras possessing finite 
dimensional representation space. The linear entropy that measures the degree of entanglement at 
the output of the beam splitter is analytically obtained. We find that the entanglement is not only 
, strongly dependent on the Hilbert space dimension but also quite related to strength the param- 

\^ I eter ensuring the temporal stability of the phase states. Finally, we discuss the evolution of the 

On ' entangled phase states. 
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1 Introduction 



Probably the entanglement phenomenon, often named quantum non-locality pji2j, contains one of 
the most interesting features of quantum mechanics. Actually, it is at the heart of current develop- 
ment of quantum information processing such as quantum teleportation [3] , superdense coding [1] and 
telecloning [5j. Entanglement also plays a key role in secure communication, for instance the Ekert 
protocol [6j based on entangled states is more robust than BB84 one [7]. It known that in quantum 
computation, the qubits are massively entangled. The preparation and characterization of entangled 
optical as well as atomic states has been studied extensively. In this respect and in a recent experi- 
mental advance, the polarization-entangled photons were generated using type I or type II parametric 
down conversion 0. 

The production of entangled states belonging to an infinite dimensional Hilbert space can be 
achieved also by adopting the standard technique of parametric down conversion [9] . The beam splitter 
is also one of the few experimentally accessible devices, which may be used to generate entangled 
states. In this sense, different developments have been reported on the analysis of beam splitter as 
entangler [lOl [TTl [121 ttS O US]. In particular, the effect of beam splitting on the spin (or SU(2)) 
coherent states for a single mode field was investigated in [14j. In the same spirit the entanglement 
via a beam splitter of SU{1,1) coherent states (Barut-Giradello |16] and Perelomov ones [T^) was 
discussed in [15]. The investigation of entanglement properties of coherent states [181 [191 [201 [21] is 
mainly motivated by the fact that the entangled nonorthogonal states also play an important role in the 
quantum cryptography [22j and quantum information processing [23]. Experimentally, the quantum 
optical systems are extensively investigated in order to generate, characterize and understand the 
entanglement properties. 

In this work, we focus on the analysis of the entangled phase states of single modes of the elec- 
tromagnetic field. It is well-known that the usual way to quantize these single modes is through the 
harmonic oscillator techniques with an infinity of states. In 1989, Pegg and Barnett [M| suggested to 
truncate up to some finite, but arbitrarily large, order the infinite dimensional representation space of 
the harmonic oscillator algebra. This was done to get rid of the difficulty related to the infinite dimen- 
sional character of the representation space, which constitutes a drawback in defining a phase operator 
in a consistent way [25ll26t [^. Here we introduce a generalized version of the Weyl-Heisenberg algebra 
that allows us to achieve our goal. We particularly consider one algebra possessing finite dimensional 
bosonic representations. This is essential in defining the Hermitian phase operator and corresponding 
temporally stable phase states. For this purpose, we use a technique based on an approach developed 
in |28j . Furthermore, we show that the phase parameter ensuring the temporal stability of phase 
states plays a crucial role in the present analysis. Finally, we deal with the entanglement of phase 
states when passed trough a beam splitter. 

The outline of the paper is as follows. In section 2, we introduce a generalized Weyl-Heisenberg 
algebra, which extends the dynamical symmetry of the usual harmonic oscillator. Subsequently, we 
discuss the corresponding finite dimensional Hilbertian representation and derive the temporally stable 
phase states, which are obtained to be dependent of ip called phase parameter. This latter allows us 
to take into account the nature of the spectrum of the system. Because of the absence of ip, one 
can immediately notice that the SU{2) phase states derived in [2U] are identical to those obtained by 
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Pegg and Barnett [21] for truncated harmonic oscillator. Furthermore, we analyze the basic features 
of the phase states. In section 3, we examine the entanglement resulting from the action of a beam 
splitter on the phase states. To investigate the degree of bipartite entanglement of the phase states, 
we determine the linear entropy. Finally, we close by some concluding remarks. 

2 Finite Fock space for generalized Weyl— Heisenberg algebra 

A basic ingredient that will be used in the forthcoming analysis is the generalized Weyl-Heisenberg 
algebra. This is generated algebraically by three elements denoted by {a^ , , N} , which are satisfying 
the commutation relations 

[N,a-]=-a-, [iV,a+]= + a+, [a", a+] = G(iV) (1) 

where G{N) = [G{N)]^ is a Hermitian function of the number operator N. Clearly, by requiring that 
G{N) = I, where I is the unity operator, we end up with the usual harmonic oscillator algebra. Also, 
for G{N) = aN + b, with two real parameters a and b, ([1]) reproduces the Wk algebra discussed in the 
context of fractional super symmetric quantum mechanics [30] . It is also important to stress that this 
algebra covers the extended harmonic oscillator worked out in |3H [52] . 

Let us consider the abstract Fock representation of the above algebra through a complete set 
of orthonormal states {|n),n G N} those are eigenstates of the number operator N, N\n) = n\n). 
In this representation, the vacuum state defined as a^\0) = and the ortho-normalized eigenstates 
constructed by successive applications of the creation operator a+. Indeed, we define the actions of 
creation and annihilation operators as 

a-\n) = /FV)e*[^(")-^("~i)]^|n- 1) 

a+|n) = /FVTl)e-^[-^("+i)-^(")]^|n + l) (2) 

where the structure function F{.) is an analytic function, with the properties F{0) = and F{n) > 
for 71 = 1, 2, • • •. The phase parameter ip will be discussed in the next by emphasizing its play a crucial 
role in constructing the phase states. In what follows we shall denote the Fock space as J- where the 
operators a"*" and a~ are mutually adjoint, a'^ = {a")^ on T. It is easy to check that F{.) satisfies the 
recursion relation 

F{n + l)-F{n) = G{n) (3) 
which gives by simple iteration the form 

n-l 

F{n) = G{m). (4) 

»n=0 

In the forthcoming analysis, we restrict ourselves to generalized oscillator algebra defined through 
the structure functions those fulfilling the condition 

F(2s + 1) = (5) 

where 2s is a positive integer value. It this case, it follows that the creation and annihilation oper- 
ators satisfy the nilpotency relations (a~)^**^^ = (a"^)^*"*"^ = 0. This means that the corresponding 
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representation is (2s + l)-dimensional. It is interesting to note that by using one can write the 
condition ([5]) as 

Tr G = 0. (6) 

where the trace is over the (2s + l)-dimensional Fock space. This new algebra is covering the following 
results: 

• (i)- The truncated harmonic oscillator introduced by Pegg-Barnett [24j : 

F{N) = N, G(Af) = 1- (2s + l)|2s)(2s|. (7) 

• (ii)- The finite dimensional oscillator algebra (k = — l/2s < 0) defined in [28| : 

F{N) = N[l + k{N -1)], G{N)=I + 2kN. (8) 

• (iii)- The truncated generalized oscillator algebra Ak {k > 0) [28] : 

F{N) = N[1 + k{N -I)], G{N) = I + 2kN - F{2s + l)\2s){2s\. (9) 



The eigenvalues of the operators G{N), corresponding to the special cases (i), (ii) and (iii) are respec- 
tively given by G{n) = 1 — (2s + l)Sn,2s, G'(n) = 1 + 2Kn and G{n) = 1 + 2«;re — F(2s + l)6n,2s for 
n = 0, 1, • • • , 2s. At this level, we have different comments in order. It is easily seen that for «; = 0, 
the case (iii) reduces to (i). As mentioned before, the algebra ([T]) is more general and covers many 
other variants of generalized harmonic oscillators. The particular cases (i), (ii) and (iii) are interesting 
in the context of quantum optics. Indeed, the case (i) was introduced by Pegg and Barnett [23] to 
define, in a consistent way, the phase operator. The others cases (ii) and (iii) constitute an extension 
of the ideas developed in [24] to define phase operator for systems with nonlinear spectrum [28| . 

Using the algebraic structure of the generalized oscillator algebra ([1]), one can introduce an operator 
that generalizes the Hamiltonian a~^a~ for the one-dimensional harmonic oscillator. Starting from 

a+a-\n) = F{n)\n) (10) 

it is easy to realize the required Hamiltonian as 

H{N) = F{N) = a+a-. (11) 

where H{N) can be regarded as the Hamiltonian corresponding to the algebra ([1]). As usual, one can 
proceed with the eigenvalue equation 

H{N)\n) = F{n)\n) (12) 

to get the solutions of the energy spectrum of a quantum dynamical system described by F{N). As 
illustration of the realized Hamiltonian, we give the explicit forms for different cases mentioned above. 
Indeed, we have for case (i) 

2s 

H{N) = ^n\n){n\ (13) 

n=0 
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and case (ii) 

2s 

H{N) = —{2s + l-n)\n){n\ (14) 

as well as case (iii) 

2s 

H{N) = ^n[l + K{n-l)]\n){n\. (15) 

n=0 

3 Temporally phase states 

We show how to build the phase states, which are temporally stable under time evolution. For this, 
recall that our Hilbert space has (2s + l)-dimensions where the actions of the operators a~ and 
on are given in ([2]). These are supplemented by the condition 

a+|2s)=0 (16) 

which can easily be deduced from the calculation of expectation value {2s\a~a~^\2s). 

As F{N) is a positive definite operator on the finite dimensional Fock space, let us consider the 
decomposition of the annihilation a~ and creation a"*" operators given by [28] 



a- =E ^/F{N), a+ = ^/F{N) (17) 
and one can show that the operator E satisfies 

£;|n) = e*!^^")-^^"-^)]'^!™- 1), n = 1,2, •••,2s (18) 
which is valid modulo 2s + 1. For n = 0, we have 

E\^) = ei[^(0)-F(2.)]^|2s). (19) 

It follows that for E\ one can obtain 

E^\n) = e-i[-f'("+i)-^Wl^|n + 1) (20) 

where n + 1 should be understood modulo 2s + 1. Note that from above relations, one can easily check 
that E is an unitary operator. Therefore, it is clear now that (jl7p constitutes a polar decomposition 
of the operators a~ and . 

To deduce the phase states associated with the finite dimensional algebra defined in the previous 
subsection, let us in the beginning look at the eigenstates of the operator E. In doing so, we solve the 
eigenvalue equation 

2s 

E\z) = z\z), \z) = Cnz''\n) (21) 

n=0 

where z G C. According to the method developed in [28], it is easy to see that the eigenvalues z should 
satisfy the discretization condition 



^2s+l 



1 (22) 
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and therefore the complex variable z is a root of unity, such as 



z = q'^ m = 0,1, •••,2s (23) 
where the q parameter is defined as 

g:=e2-V(2^+i). (24) 
Thus, these tell us that the normalized eigenstates \z) = \m,ip) oi E are of the form 

|m, 99) = , y e-^-^(")^g'""|n) (25) 

* n=0 

where the parameters m G Z/(2s + 1)Z and 99 G M. It is easy to see that 

E\m,^) =e'^^\m,^), 6^ = m—— (26) 

which reflects that is indeed a phase operator and therefore \m,ip) are the required phase states. 
Note that, in the particular case 99 = 0, the states |m,0) correspond to an ordinary discrete Fourier 
transform of the basis {|n) : n = 0, 1, • • • , 2s} of the (2s + l)-dimensional space J-. We show that the 
phase states \m,{p) are temporally stable under "time evolution", namely 



e 



-'^'^^^'Im^if) = \m,ip + t) (27) 



for any value of the real parameter t. The parameter ip plays a key role to ensure the temporal stability 
of the states \m,ip). For fixed ip, they satisfy the equiprobability relation 

\{n\m,(f)\ = , ^ n,mGZ/(2s + l)Z. (28) 
V2s + 1 

They constitute an orthonormal set 

{m,if\m' ,ip) = 6m,m' m, m' G Z/(2s + 1)Z (29) 
for fixed ip and satisfy the closure property 

2s 

^ \m,ip){m,ip\ = I. (30) 

m=0 

Finally, it is interesting to note that the temporally stable phase states are not all orthogonal. Indeed, 
the overlap between two phase states \m',ip') and \m,ip) reads as 

1 

{m, ip\m', = (31) 

n=0 

where the function p is 

2s + 1 

p{m — m', (p — 99', n) = —{m — m')n H (99 — ip')F{n) (32) 

and q is given in 
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At this level, we would like to emphasize the key role of the phase parameter 99 introduced starting 
from relations ([2]) that defines the actions of creation and annihilation operators. Indeed, firstly 
ip ensures the temporal stability of the phase states \m,ip) (j25p under time evolution. Secondly, if 
one ignores if, i.e. (/9 = 0, \m,ip) reduce to those derived by Begg and Barnett using the truncated 
harmonic oscillator [21] . This means that for two generalized Weyl-Heisenberg algebras characterized 
by different structure functions F{N), we end up with the same phase states. Thus, to differentiate 
between these states, the presence of if is necessary in their forms. In this sense, it must be noticed 
that the SU{2) phase states obtained in [29] are similar to those of Pegg and Barnett despite the fact 
that the involved symmetries are different. Hence, to avoid such problem we introduced ip that allows 
us to keep the trace of the symmetry and dynamics of the system under consideration. Note also that 
(p plays a key role in relating phase states and mutually unbiased bases [28j. 



4 Beam splitting of phase states and entanglement 



We briefly describe the effect of beam splitting on a known input state with the vacuum at the second 
port. We assume that the vertical input beam is always prepared in the Fock ground state and 
the state of interest is the input state in the horizontal beam, which by convention precedes that of 
vertical beam. Algebraically, the beam splitting transformation is described as follows. The input field 
described by the usual harmonic oscillator operator is superposed on the other input field with 



operator 05 by a lossless symmetric beam splitter, with the transmission t and reflection r coefficients. 

and 



The output field operators and are given by 



a 



± 



B{o)4B\e) 



where the unitary beam splitter operator of angle 9 is 

B{e) = exp 

The action of B{9) on the state |n,0) reads as 

71 f 

B{e)\n,Q) = Y, 



-9 {afa^ + a^) 



(33) 



(34) 



n! 



p=0 



p\{n — p)l 



tP{irY~P\p,n-p) 



(35) 



where the transmission and reflection coefficients are defined by 



t 



cos ■ 



sm ■ 



(36) 



2 2 

It is well known that among pure states of a single mode field, only harmonic oscillator or Glauber 
coherent states do not become entangled upon beam splitting. Indeed, for the 50:50 beam splitter 
with an ordinary Glauber coherent state incident on one port and a vacuum on the other, the beam 
splitting result is a product of two Glauber states. This procedure do not provides us with an entangled 
system. In general, any other pure state at the input results in an entangled state at the output. This 
explains how the beam splitter acts as entangler and why it is used in experiments as device to generate 
entangled states. 
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Now, let us proceed in our case to rewrite the action of the beam sphtter operator on the phase 
states ([25]) as 

^ 2s n j j 

S(^)|m,,.)^|0) =-^_ j;^;^^^^^ g-t-(.rr- e-^^W- \p,n-p). (37) 

The corresponding output density is then given by 

pi,2 = B{9)\m, if) |0)(m, ip\ {0\B\9) (38) 

It fohows that one can write the bipartite reduced density matrix pi2, which is obtained by tracing 
out the second system. This is 

pi = Tr2/9i,2 (39) 

where label 2 is reflecting the trace on second states. It is easy to see that the reduced density matrix 
is given by 

2s 2s mm{2s-n,2s—n') 

P^ = J2Y1 E c(n,0^(<T)|n)(n'| (40) 

n=On'=0 /=0 

where the coefficients c{n, I) have the form 

c(n, /) = 1 \ (?™("+')r(ir)' e-*^("+')^. (41) 

Next, we are interested in the amount of entanglement of the beam splitter output states. In doing 
so, we examine the entanglement of phase states when passed through a beam splitter by using the 
linear entropy [33] 

S = l-Tr(p?). (42) 
Clearly, it goes to zero for a pure state. After some algebra, we show that (|42p rewrites 

2s 2s min(2s— n,2s— n') min(2s— n,2s— n') 

^=1-I;E E E s{n,n',l,l') (43) 

n=On'=0 1=0 l'=0 

where s{n, n', I, I') is given by 

1 , f2{n+n') 2(1+1') 

s{n,n',l,l') = + ^)2 ^"^'^^"'" '''' V(n + mn' + I'V-in + mW + 0'^^ (44) 

and the phase term reads as 

^{n, n', /, /') = [F{n + I) + F{n' + I') - F{n' + I) - F{n + /')] (45) 

We have some remarks in order. Firstly, according to last equations we notice that the linear entropy 
is m-independent. Secondly, one can verify that the functions (/>(n, n', /, ?') satisfy some symmetries 
with respect to interchange of their quantum numbers. These are 

(/)(n, n', /,/') = —(/)(n, n', /), (p{n,n' = (p{n' ,n,l,l') (46) 

which can be used to express the linear entropy in terms of cos (p{n,n' ,1,1') instead of e^*'^("'"'''''') 
together with the symmetry of the summations in (|43|) . 
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For the temporally stable phase states considered in section 3, we study numerically the behavior 
of the first-order entropy as a function of the relevant parameters s, (p and the reflection coefficient r 
as well. We shall focus on the phase states associated with the generalized Weyl-Heisenberg algebra 
defined through the structure function 

F(iV) = ^(2s + l-iV). (47) 

We first plot the entanglement against (p and for different finite dimensional Hilbert of (2s + 1)- 
dimensions. 




s 



Figure 1: Linear entropy of phase states as a function of i? = and 99 for s = 1/2. 
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Figure 2: Linear entropy of phase states as a function of i? = and (/? for s = 1. 




According to the above figures, we note that for all value taking by s there is a maximum entanglement 
at the point = 1/2, which corresponds to a 50:50 beam splitter. Furthermore, we have two 
interesting features in order. The first is that for qubits s = 1/2, the degree of entanglement is (p- 
independent as shown in Figure 1. However, for qutrits (s = 1) and fixed reflection coefficient, the 
behavior of the linear entropy is Gaussian. The maximum is reached for = 1/2 and ip = n, see 
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Figure 2. This changes completely for higher dimensional Hilbert spaces. Indeed, for s = 3/2, for 
fixed value of r^, the linear entropy undergoes an initial rapid increase to reach the maximum entropy 
for if = TT followed by a slower drop and then increases when ip does go to 27r as it shown in Figure 3. 

As we explained above, the parameter (p plays a key role in defining the phase states and ensuring 
the temporal stability of phase states. It follows that it is interesting to examine the linear entropy in 
term of ip for s fixed and = 1/2 for which one has the maximum of entanglement of phase states. 




Figure 4: Linear entropy of phase states after being passed through a 50:50 beam splitter against if for 

5 = 1/2,1, •••,5. 

The behavior of linear entropy as function of the parameter ip is represented in Figure 4. It is easily 
seen that for s = 1/2, the linear entropy is constant. For qutrits, i.e. s = 1, the entropy looks like 
a gaussian with a maximum around ip = 2tt. For higher levels quantum systems, the linear entropy 
behaves differently. For instance, with phase states of quartits, i.e. s = 3/2, the degree of entanglement 
increases for ip < it, decreases for it < ip < 3tt/2 and increases after. As the parameter (p is deeply 
related to time evolution of phases states, the figure 4 can be viewed as representing the temporal 
evolution of the degree of entanglement of phase states. 

In Figure 5 we plot the linear entropy of phase states versus the parameter s for different values 
of p to underline its basic behavior. 



3 TT 



20 
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Figure 5: Linear entropy of phase states after being passed through a 50:50 beam splitter against s for 

different values of (p. 

Following this figure, we notice that the linear entropy initially rapidly rises for lower s but increases 
slowly for s > 10. However for higher s, the linear entropy does not approaches zero. This shows 
that the degree of entanglement is increasing when the dimension of the generalized Weyl-Heisenberg 
becomes large. 

5 Concluding remarks 

The main idea of the present work is the investigation of the degree of entanglement of temporally 
phase states. The latter are constructed as cigcnstatcs of a unitary phase operator resulting from 
the polar decomposition of ladder operators of finite dimensional Weyl-Heisenberg algebra. Using the 
beam splitter as entangler, we computed the degree of entanglement of phase states by mean of linear 
entropy. We investigated the behavior of this quantity as function of the dimension of the Hilbert 
space, the transmission-reflection coefficients of the beam splitter and the parameter ensuring the 
temporal stability of the phase states. 

It is clearly shown that the maximal entanglement is provided by 50:50 beam splitter. Also, the 
entanglement increases with increasing dimension of the system. For instance the phase states for 
qutrits are more entangled than qubits. The entanglement of phase states when passed through a 
beam splitter is strongly dependent on the parameter ip which play an essential role in defining the 
eigenstatcs of the unitary phase operator and ensure their temporal stability. 

Finally, we note that we did not consider the limiting case s — > oo. In fact for infinite dimensional 
Hilbert space, one can not define unitary phase operator in a consistent way. Consequently the phase 
states associated with the ordinary Weyl-Heisenberg can not be constructed using the formalism 
developed in this work. 
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